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Abstract. Magnetoelastic properties of the spin-1/2 Ising-Heisenberg model on doubly decorated planar 
lattices partially amenable to lattice vibrations are examined within the framework of the harmonic ap- 
proximation and decoration-iteration transformation. It is shown that the magnetoelastic coupling may 
lead to a spontaneous distortion of the vibrating decorating atoms and the mutual interplay between quan- 
tum spin fluctuations and local lattice deformations enhances typical quantum features like the quantum 
reduction of the magnetization in the ground state of the quantum antiferromagnetic phase, while it does 
not affect the ground-state behaviour of the classical ferromagnetic phase. It also turns out that the spon- 
taneous distortion is responsible for a much more pronounced reduction of the critical temperature in the 
quantum antiferromagnetic phase than in the classical ferromagnetic phase. 

PACS. 05.50.-|-q Lattice theory and statistics - 05.70.Jk Critical point phenomena - 75.10.-b General 
theory and models of magnetic ordering - 75.30.Kz Magnetic phase boundaries - 75.40.Cx Static properties 
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1 Introduction 

Classical and quantum spin models traditionally attract 
an appreciable scientific interest as they bring an insight 
into various aspects of phase transitions, critical phenom- 
ena and many other diverse phenomena of a cooperative 
nature |1I2I3I4| . However, the predominant number of sug- 
gested lattice-statistical spin models tacitly assumes that 
atomic spins are placed at rigid lattice positions and they 
interact among themselves by means of a constant (usually 
nearest-neighbour) exchange interaction. This conjecture 
represents an obvious idealization, because atoms are car- 
rying out persistent oscillatory motion around their equi- 
librium lattice positions even at zero temperature and the 
exchange interaction may strongly depend on the inter- 
atomic distance ^ . Early attempts to deal with magnetoe- 
lastic properties of compressible spin models were based 
on a rather simple but unrealistic assumption that the ex- 
change interaction between the atomic spins depends just 
on the average volume |6I7I8| . To take into account the 
effect of fluctuations of lattice spacings on a magnitude of 
the exchange interaction, the more realistic compressible 
spin models were later proposed on assumption that the 
exchange interaction varies linearly with an instantaneous 
interatomic distance .!) II) I I 12 13, 14) . Among other mat- 
ters, it has been shown that the spin-lattice interaction 
may change a continuous (second-order) phase transition 
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to a discontinuous (first-order) one [15116) , provoke an in- 
teresting renormalization of critical exponents |14ll7j . or 
cause the lattice instability with respect to a spontaneous 
spin-Peierls dimerization jlS 19 20 21 22 . 

It is worthy to notice that almost all exactly tractable 
spin models including the spin-lattice interaction are of a 
classical nature such as the Ising models (9il0llll23l24| or 
the classical Heisenberg models [25126127128129] . Despite 
a considerable effort, rigorous treatment of quantum spin 
models with the spin-lattice interaction is usually incred- 
ible task except of a few highly valuable cases like the 
spin-1/2 XX chain in a transverse magnetic field display- 
ing the spin-Peierls phenomenon [TS 19 20 21 22 or the 
compressible version of the spin-1/2 XX chain in a trans- 
verse magnetic field |30j. Accordingly, one usually has to 
rely on application of some approximative method when 
treating a quantum spin model with the spin-lattice inter- 
action [31 32 33 34 35 36 37 38 . 

An important challenge for theoretical studies currently 
represents an investigation of the effect of the spin-lattice 
interaction on ground state, phase transitions and criti- 
cal phenomena of quantum spin models. Bearing this in 
mind, it is also highly desirable to search for simpler ex- 
actly tractable quantum spin models with the spin-lattice 
interaction, which could provide a deeper understanding 
of the effect of a mutual interplay between quantum spin 
fluctuations and fluctuations of lattice spacings on a crit- 
ical behaviour. The hybrid Ising-Heisenberg models de- 
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fined on planar lattices, which are occupied in part by the 
classical Ising spins and partly by the quantum Heisen- 
berg spins, belong to a few eligible candidates capable 
of this study [3914014 1 42 43 44 45 . It is noteworthy that 
the exactly solved Ising-Heisenberg planar models have 
already proved their usefulness in predicting several un- 
usual ordered or disordered quantum phases [39II40| 41I|42I , 
explaining a partial lifting of the ground-state degeneracy 
in disordered spin-liquid phases |41I42| , providing a deeper 
insight into a weak-universal critical behaviour [43144] or 
quantum correlations from the viewpoint of quantum com- 
putation science [IS]. The main objective of the present 
work is to extend the spin-1/2 Ising-Heisenberg model 
on doubly decorated planar lattices from Ref. [31] by ac- 
counting for the spin-lattice interaction treated within the 
harmonic approximation. Our calculation will closely fol- 
low the approach invented by Bellucci and Ohanyan [46) 
for the spin-1/2 Ising-Heisenberg sawtooth chain with the 
spin-lattice interaction except that the spin degrees of free- 
dom of the Heisenberg spins will be treated before inte- 
grating out the lattice degrees of freedomQ 

Finally, it is worthwhile to remark that the magnetic 
structure of the model under investigation has been in- 
spired by two layered metal frameworks of polymeric coor- 
dination compounds 2Rh2{CH3COO)4 • K3Co(CN)6 [17] 
and [Cu2(L)]3[Co(CN)6]2 • 6CH3OH • 2DMF gg in which 
metal ions are arranged so as to build the doubly deco- 
rated square lattice and the doubly decorated honeycomb 
lattice, respectively. Unfortunately, the trivalent Co"^"'' ions 
are diamagnetic in both aforementioned compounds due 
to a strong ligand field of the cyano group and hence, there 
does not exist nowadays any appropriate results of mag- 
netic measurements to compare with our theoretical pre- 
dictions. However, the recent progress in a targeted design 
of molecular-based magnetic materials with desired crys- 
tal structure and magnetic properties gives us hope for a 
preparation of experimental representative(s) of the model 
under investigation [49150151] . As a matter of fact, the 
rational synthesis of isostructural analogs of both afore- 
mentioned polymeric complexes 47 48, that would con- 
tain within the complex anion [M(CN)6]'^~ the low-spin 
Fe^"*" ion (5=1/2) instead of the diamagnetic Co'^"'" ion 
(5=0), represents one of the most promissing routes for a 
design of suitable model compound and might be regarded 
as a challenge for molecular engineering. Besides, it is also 
our hope that the results presented in the following may 
bring insight into many important features arising from 
the spin-lattice interaction, which could manifest them- 
selves in a wider class of insulating magnetic materials. 

The outline of this paper is as follows. In Section [51 
we first provide a rather detailed description of the model 
under investigation together with the most crucial steps 



^ Note that the reverse procedure as used in Ref. 46 is not 
actually exact, because one cannot simply factorize the bond 
partition function into two parts (one of which is independent 
of the distortion parameter, while the other one does depend 
on the distortion parameter) before integrating out the lattice 
degrees of freedom due to a non-commutability of both parts 
(perturbed and unperturbed) of the same bond Hamiltonian. 



of the calculation procedure, which enables us to explore 
ground-state and finite-temperature phase diagrams, the 
relevant behaviour of order parameters and the distortion 
parameter. The most interesting results are then presented 
and detailed discussed in Section [HJ Finally, the summary 
of the most important scientific achievements is mentioned 
with several concluding remarks in Section [4] 

2 Model and its solution 

Consider the spin-1/2 Ising-Heisenberg model defined on 
doubly decorated planar lattices as schematically illus- 
trated in Fig. [T] on a particular example of the doubly 
decorated square lattice. Within the framework of this 
model, each nodal site of some doubly decorated planar 
lattice is occupied by the Ising atom with the spin a = 1/2 
while each its bond incorporates a couple of the quantum 
Heisenberg atoms with the spin 5 = 1/2. Suppose further- 
more that the Ising atoms are placed at rigid lattice po- 
sitions in contrast with the decorating Heisenberg atoms, 
which are relaxed from the condition of a perfect lattice 
rigidity and may consequently change their equilibrium 
lattice positions after taking into account the elastic en- 
ergy penalty for each distorted bond within the harmonic 
approximation. It is noteworthy that the assumption of 
lattice rigidity of the Ising atoms perhaps represents the 
most important drawback of our model, which is however 
quite reasonable if: (i) the mass of the Ising atoms is much 
larger than the mass of the Heisenberg atoms; (ii) the 
model is defined on doubly decorated planar lattices with 
a high enough coordination number of the Ising atoms (the 
displacement of the Ising atoms from their equilibrium lat- 
tice positions would then cost much larger amount of the 
elastic energy owing to a deformation of greater number 
of bonds); (iii) the spring stiffness constant of the Ising 
bonds is much stronger than the elastic constant of the 
Heisenberg bonds; (iv) various combinations of previous 
three particular cases. 

For further convenience, let us write the total Hamil- 
tonian of the investigated model system as a sum of bond 
Hamiltonians 

Nq/2 Nq/2 

^= E E(^r+^f)- (1) 

i=l 1=1 

Each bond Hamiltonian Hi, which is further splitted into 
the magnetoelastic part and the pure elastic part , 
contains all the interaction terms of one couple of the 
Heisenberg spins from the zth bond (see Fig. [1]) [N denotes 
the total number of the Ising spins and q determines the 
number of their nearest neighbours]. The magnetoelastic 
part of the bond Hamiltonian is given by 

- - [J + A(xa + x,2)][Zi(5ri5f2 + S^SD + S^^S^^ 

- (Ji - AiX,i)5fi(Tfi - (Ji - AiX,2)5f2fTf2. (2) 

Here, Sj^ and a^^ (a = 1,2, 7 = x,y,z) denote spatial 
components of the standard spin-1/2 operator, the pa- 
rameter J labels XXZ interaction between the nearest- 
neighbour Heisenberg spins, Z\ is a spatial anisotropy in 
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# Ising spin ct=1/2 

# Heisenberg spin 5= 1/2 

Fig. 1. The segment from the doubly decorated square 
lattice. Blue circles label rigid lattice positions of the Ising 
spins a = 1/2, while the red ones denote lattice positions 
of the decorating Heisenberg spins 5 = 1/2 prone to lattice 
vibrations. 

this interaction and the parameter Ji stands for the Ising 
interaction between the nearest-neighbour Ising and Heise- 
nberg spins. Furthermore, the parameter Xia denotes a 
displacement of the Heisenberg spin Sia from its equi- 
librium position referred relative to its nearest-neighbor 
Ising spin CTm, whereas Xia > {xia < 0) corresponds 
to an elongation (contraction) of the relevant distance. It 
is worthwhile to remark that only a rather small lattice 
displacements are feasible within the harmonic approxi- 
mation, which consequently allows one to perform a lin- 
ear expansion of the distant-dependent exchange interac- 
tions around their non-distorted values J and Ji. Hence, 
it follows that a linear increase (decrease) in the Heisen- 
berg and Ising interaction originates from a contraction 
(elongation) of a relevant bond length and this distortion- 
induced change in the exchange interactions J and Ji is 
proportional to the magnetoelastic (spin-lattice) couplings 
A and Ai, respectively. Finally, the pure elastic part of the 
bond Hamiltonian takes into account a kinetic energy of 
both Heisenberg atoms prone to lattice oscillations as well 
as the elastic energy penalty for one Heisenberg bond and 
two Ising bonds to be treated as classical harmonic oscil- 
lators 

nt = + f (.a + x.,f + ^ (4 + ^h) . (3) 

Above, pia — denotes the momentum of the Heisen- 

berg atoms with the mass M, K and Ki label spring stiff- 
ness (bare elastic constants) of the Heisenberg and Ising 
bond, respectively. The first term thus represents the ki- 
netic energy of both the Heisenberg atoms, the second 
term stands for the elastic energy of the Heisenberg bond 
and the last third term represents the elastic energy of 
both Ising bonds. Last but not least, it is worthy to note 
that the elastic part of the bond Hamiltonian Q takes into 



consideration only 'longitudinal' lattice deformations in 
directions of lattice bonds, which should provide the most 
important effect because in-plane or out-of-plane 'trans- 
verse' lattice deformations would couple quadratically in 
the Hamiltonian and thus, these less important higher- 
order terms would just effectively cause rescalling of the 
elastic term only. 

The crucial step of our procedure represents calcula- 
tion of the partition function. Different bond Hamiltonians 
obviously commute with each other, i.e. [HijUj] = 0, and 
in the consequence of that, the total partition function can 
be calculated when summing the product of bond parti- 
tion functions over all possible spin configurations of the 
Ising spins 

Nq/2 

^=En^- (4) 

Each bond partition function Zi in fact represent the local 
partition function of one couple of the Heisenberg spins 
from the ith bond and can be defined as follows 

oo oo 

^JJOx^ JJOx^ exp(-/37^f )Tr, exp(-/3^r). 

— OO — OO 

(5) 

Here, (3 = 1/(^8^), is Boltzmann's constant, T is the 
absolute temperature, h is Planck's constant and the sym- 
bol Tri means a trace over spin degrees of freedom of the 
ith couple of the Heisenberg spins. 

To proceed further with calculations, it is advisable to 
perform local canonical coordinate transformation unam- 
biguously given by 

Xii = (y»i + 2/12) , ^ ^^'^ ~ ^''^^ ' 

or its inverse 

Vil = ^{Xil + Xt2) , yi2 = - X^2)] (7) 

which introduces two new spatial coordinates yn and yi2 
instead of local displacement parameters xn and Xi2. The 
spatial coordinate yn is connected to an average displace- 
ment of both Heisenberg spins from their equilibrium lat- 
tice positions, while the other spatial coordinate yi2 re- 
flects a possible assymetry in the displacement of both 
Heisenberg spins relative with respect to their nearest 
Ising neighbors. The main advantage of the canonical co- 
ordinate transformation ^ lies in a separable form of the 
potential energy from the elastic part of the bond Hamil- 
tonian (which subsequently allows an independent inte- 
gration over both new spatial coordinates yn and 2/^2) 
while keeping the kinetic term unchanged. Applying the 
canonical coordinate transformation (|6]), the elastic part of 
the bond Hamiltonian ([3]) can be actually rewritten using 
the new coordinates yia and their conjugated momenta 
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M-^ into the following form 



2M 



(8) 



A substitution of the elastic part of the bond Hamiltonian 
([5]) into Eq. ((S]) yields the following result for the bond 
partition function after a straightforward integration over 
both momenta 



2ttM 



I3K 



dyn exp 



X /d2/,2exp ( -^K^yj^ Tr, exp(-/?Hr). (9) 



Tr,exp(-/3Hr) = exp 



which are explicitly given by the following formulas 



Eii.. 



(10) 



(11) 



i?z3,,4 = f ±^V^U^^^J2^^)V"(^^. (12) 

For the sake of brevity, the four energy eigenvalues pT|) - 
are expressed in terms of the distant-dependent ex- 
change interactions 

J = J + A(a;,i +xa) = J + V2Ay,i, 



Ji = Ji- AiXii = Ji - (ya + yi2) , 
J2 = Ji- AiXi2 = Ji - ^ (ya - yi2) ■ 



(13) 



It should be noted here that the former couple of eigenval- 
ues (|lip varies linearly with the displacement parameters 
yn and yi2 in contrast to the latter couple of eigenvalues 
involving both these displacements inside a square 
root, which would consequently preclude an integration 
over lattice degrees of freedom at the level of the bond par- 
tition function. However, the latter two eigenvalues (|12l) 
can still be expanded into a power series in terms of both 
displacement parameters yn and yi2, whereas it is suffi- 
cient to retain only the terms up to the first order and 
neglect all higher-order terms on behalf of small magni- 
tude of both displacements. Using this approach, one eas- 
ily obtains four eigenvalues of the Hamiltonian "Hf* that 



linearly depend on the displacement parameters yn and 
yi2 according to the relations 



Eii.. 



-4Ty(a,i+a,i)±^Ai(a,, 

yki 



[At^i«i+<i)], 



2^/2 



2 



It is quite evident from Eq. ^ that it is now necessary 
to find the partial trace over spin degrees of freedom of 
the ith couple of the Heisenberg spins in order to get ex- 
pression suitable for a subsequent integration over spatial 
coordinates. This can be rather easily accomplished by di- 
agonalizing the magnetoelastic part of the bond Hamilto- 
nian ^ and considering a trace invariance, which enables 
one to express this partial trace in terms of four eigenval- 
ues of the Hamiltonian 'H™ 



yki 
2V2 



A± 



JiAi{cjf^ + af^f - 2J AA^ 



{JAf 



(14) 



(15) 



At this stage, the four energy eigenvalues given by Eqs. ([T4f 
and can be used for calculating the partial trace in 
the magnetoelastic part of the bond partition function 
(|10p and the resultant expression can be further substi- 
tuted to the bond partition function in order to per- 
form an integration over lattice degrees of freedom. After 
a straightforward but little bit cumbersome integration 
over the spatial coordinates yi\ and yi2 one finally gets 
the bond partition function 



2 exp 



16(^1+2^) 



P^W-ijj\u:2 
(3Al{at,-at^f 



exp 



4 



exp 



4 



cosh — ( 4 Ji 



lQ{Ki + 2K) 
AAi 



Ki + 2K^ 
2 - 2JAA^]^ 



'i2 



\^{Ki^2K){Jl{al^-al^Y^[JAy 



X cosh 



/3A 



^f2) + {J^r 

JMcjl^-al^f-2JAA^ 



8{Ki + 2K) 



(16) 



whereas = Ki /M and 0J2 — a/ {Ki -f 2K)/M denote 
two angular frequencies that correspond to out-of-phase 
and in-phase normal mode vibrations of both Heisenberg 
spins. It is of fundamental importance that the bond parti- 
tion function does not already depend neither on the 
Heisenberg spin degrees of freedom nor the lattice degrees 
of freedom, but only upon the two Ising spins an and 17^2 ■ 
Accordingly, one may exploit the generalized decoration- 
iteration transformation |52l53l54l55j in order to substi- 
tute the bond partition function (I16p by a simpler equiva- 
lent expression depending just on two Ising spins an and 



Zi{aa,ai2) = i?oexp(/3i?icriicri2)- 



(17) 



The physical meaning of the generalized decoration-itera- 
tion transformation (|17p lies in replacing all the interac- 
tion parameters belonging to the zth bond Hamiltonian T-Li 
through the unique effective interaction i?i between the 
two nodal Ising spins an and ai2- This mapping transfor- 
mation should of course hold independently of spin states 
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of both these Ising spins and this 'self-consistency' condi- 
tion [52l53l54l55j unambiguously determines yet unknown 
mapping parameters -Ro and Ri 



Rq = 2exp 



= 2 In 



16{Ki+2K) 

V2 



(18) 
(19) 



which are expressed by means of the functions Vi and V2 



Vi = exp 



- exp 



V2 — exp 



T 

- + 
4 



16(^1 + 2K) 



cosh 



4(i^i + 2K) 



cosh 



2 



^AAi 



8{Ki + 2K) 
I3A^A 



16K, 



exp 



X cosh 



J! + (J Ay 



/3J 

T 

PAPi 



2K) 



8{Ki + 2K)_ 

and the auxiliary function Pi defined as 

^ _ JiAi- 2JAA^ 
' " ^Jf + {JAf ■ 



(20) 



(21) 



If one substitutes the decoration-iteration transforma- 
tion (IT7I) satisfying the 'self-consistency' condition (IT8|- 
(PT|) into the partition function one readily gains a 
mapping relationship between the partition function Z of 
the spin-1/2 Ising-Heisenberg model on the doubly dec- 
orated planar lattices amenable to lattice vibrations and 
respectively, the partition function Zjm of the equivalent 
spin-1/2 Ising model on a corresponding undecorated and 
rigid planar lattice 



Z(/3, J, Ji,Z\,A^i,^, Ki) 



<«/'ZiM(/3,i?i). 



(22) 



Exact expressions for the partition function of the spin- 
1/2 Ising model are well known for several planar lat- 
tices [56I57I58I59I60) and accordingly, the mapping rela- 
tion (|22p can be rather straightforwardly employed for 
calculating the partition function of the model under in- 
vestigation as well. Apart from this fact, it is quite evident 
from Eq. that the spin-1/2 Ising-Heisenberg model on 
a vibrating doubly decorated planar lattice becomes crit- 
ical if and only if its equivalent spin-1/2 Ising model on a 
corresponding undecorated (rigid) lattice becomes critical 
as well. Critical frontiers of the spin-1/2 Ising-Heisenberg 
model on a vibrating doubly decorated planar lattice can 
be thus easily obtained from a comparison of the effec- 
tive coupling dT51) of the equivalent spin-1/2 Ising model 
with its critical value. For instance, the critical lines of 
the spin-1/2 Ising-Heisenberg model on a vibrating dou- 
bly decorated square lattice shown in Fig. [1] can be deter- 
mined from the comparison of the effective coupling (jl9p 
with the relevant critical point of the spin-1/2 Ising square 
lattice [56l 



\Ri\ 



21n(l + \/2). 



Here, Tc denotes the critical temperature and the absolute 
value |i?i| reflects the same value of the critical temper- 
ature irrespective of whether the effective interaction is 
being ferromagnetic i?i > or antiferromagnetic i?i < 0. 

Last but not least, let us derive expressions determin- 
ing the behaviour of the order parameter, which might 
be either the spontaneous magnetization for a particu- 
lar case of the ferromagnetic long-range ordering or the 
spontaneous staggered magnetization for a special case of 
the antiferromagnetic long-range ordering. The uniform 
as well as staggered spontaneous magnetizations of the 
Ising spins tua and sa can readily be derived by employ- 
ing exact mapping theorems developed by Barry et al. 
|6ll62l63l64j 



EmiM(/3,i?i >0),(24) 
siM(/3,i?i < 0), (25) 



according to which the canonical ensemble average (• • • ) of 
any function involving only the Ising spins in the spin-1/2 
Ising-Heisenberg model on a vibrating doubly decorating 
lattice directly equals to the canonical ensemble average 
(■ • • )iM of the same function of the Ising spins in the equiv- 
alent spin-1/2 Ising model on a corresponding undecorated 
rigid lattice. For completeness, let us quote explicit formu- 
las for the uniform and staggered spontaneous magnetiza- 
tions of the Ising spins in the spin-1/2 Ising-Heisenberg 
model on the vibrating doubly decorated square lattice, 
which are expressed in terms of the relevant exact result 
for the spontaneous (either uniform or staggered) magne- 
tization of the spin-1/2 Ising model on the square lattice 
[55155] 



1 

mA — — 
2 



sa = — 
2 



1 - 



1 



1 - 



sinh* (^i?i/2) 
1 

sinh4(/3|i?i|/2) 



(fori?i>0), (26) 
(fori?i < 0). (27) 



On the other hand, it is worthy of notice that the cal- 
culation of the uniform and staggered spontaneous mag- 
netization of the Heisenberg spins and ss is somewhat 
more involved. For this purpose, one may exploit the exact 
generahzed Callen-Suzuki spin identity [67.68.69.7 71] 



TUB 



_ / •5'fl + <S'f; 



(28) 



'^il + '-'12 




exp{-f3Hi 



(29) 



1 f fdyadqn f fdyi2dqi2 



-Tr,, 



QZ QZ 



exp(-;3Hi) 



It should be nevertheless mentioned that both expressions 

(23) 

appearing inside canonical ensemble average in Eqs. (1281) 
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and (j29|l can be obtained following the same procedure as 
worked out previously when calculating the bond parti- 
tion function Zi. The spontaneous magnetization of the 
Heisenberg spins niB can in turn be related to the spon- 
taneous magnetization of the Ising spins uia through the 
exact formula 



rriB = mA 



(30) 



with the expression V\ given by Eq. ((20)) and the coefficient 
Qi defined as 



:exp 



16(^1 + 2K) 



sinh 



2 



l3AAi 



8{Ki + 2K] 



(31) 



After straightforward but a rather lengthly algebraic ma- 
nipulation, the spontaneous staggered magnetization of 
the Heisenberg spins sb can also be expressed in terms 
of the spontaneous staggered magnetization of the Ising 
spins sa through the relation 



SB 



SA_ 
V2 



J1Q2 



P2{AQ2+PlQ3) 



^Jl + {JAf MKi + 2K) 
which contains two newly defined functions Q2 and Q3 



(32) 



(^2 = exp 
X sinh 

Q3 = exp 
X cosh 



16(i^i + 2K) 



{JAf ^ 
PP? 



fiAPi 



8{Ki + 2K) 



lQ{Ki + 2K) 



1^ 



{JAY 



l3APi 



8(A'i + 2K) 

and the another auxiliary function P2 defined as 
Al JiPi 



P2 = 



^Jl + {JAY Jl + {JAY 



(33) 



(34) 



the mean displacements of both Heisenberg atoms from 
the same bond are equal to each other. For brevity, let us 
merely quote the final result for the mean displacement 
calculated by making use of the former identity 



{Xii) = {Xi2) 



(1-4£im)[A(F2 - 203)-PiQi] 



^2 8V2{Ki + 2K) 

(1 + AeiM)[A{Vi - 2Q5) - AiQi + 2AAQe] 
8Vi{Ki+2K) 



(35) 



The newly defined functions Q^-Qe entering the above 
formulas psp are given by 



Qi = exp 



"/3J 
- + 


(3AI 1 


cosh 


fiJi 


(3AAi ' 






lQ{Ki + 2K) 


2 


8{Ki + 2K) 






fiA^A" ■ 


cosh 


'PJA 


pA^A 




_ 4~ 


^4(i^i + 2if) 


2 


A{Ki + 2K) 


PJ 




sinh 


'(3JA 


13 A^ A 




_ 4~ 


A{Ki+2K)_ 


2 


A{Ki + 2K) 



Q5 = exp 



Qe = exp 



and the quantity eim = {<^ii<^i2)iM labels the pair corre- 
lation between the nearest-neighbour spins of the corre- 
sponding spin- 1/2 Ising model on undecorated rigid lat- 
tice. Hence, it follows that the foreknowledge of the nearest- 
neighbour spin correlation is needed in order to complete 
a calculation of the average displacement. For the sake of 
completeness, let us therefore quote also the relevant exact 
result for the nearest-neighbour pair correlation function 
of the spin- 1/2 Ising square lattice |2l59j with the effective 
nearest-neighbour coupling /3i?i 



£lM = 



1 



2 tanh' 



V 2 



1 



/C(fc) 



in which the function JC(k) denotes the complete elliptic 
integral of the first kind 



JC{k) = 



d0 



Vl - /fc2 sin^ 9 



with the modulus k = 2 smh{^)sedY {^) . 



Finally, the generalized Callen-Suzuki identity 
will be utilized for a calculation of the mean displacement 
of both Heisenberg atoms from their equilibrium lattice 
positions. Following the same procedure as described pre- 
viously by the calculation of the bond partition function 
Zi, the mean displacement of the Heisenberg atoms can 
be calculated with the help of the exact identities 



{y^ 



{yt2) = 



dyi2dqi2 




ynTvi e:iq){-l3Hi 



2/i2Tri exp(-^'Hi 



It can be easily proved that the latter displacement equals 
identically to zero (2/^2) = 0, which serves in evidence that 



3 Results and discussions 

In this section, we will describe in detail the most inter- 
esting findings obtained for the spin- 1/2 Ising-Heisenberg 
model on doubly decorated planar lattices whose decorat- 
ing atoms are prone to lattice vibrations. Before discussing 
the most interesting results, however, it is worthwhile to 
remark that the rigid limit {K/J — 7> 00) of the model 
under investigation has been thoroughly investigated by 
one of the present authors in Ref. [33] • The main focus 
of this work is therefore to highlight the most crucial dif- 
ferences between the magnetic behaviour of the spin- 1/2 
Ising-Heisenberg model on the perfectly rigid doubly dec- 
orated lattice |39j and respectively, the doubly decorated 
lattice whose decorating sites are relaxed from the condi- 
tion of a perfect rigidity. Another particular goal of the 
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JJJ= 
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1.0 




FLP 






QAP 
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[K/J]' 

Fig. 2. Ground-state phase diagram in the J/K — A plane 
for the fixed interaction constants Ji/J = 1.0 and A/ J — 
1.0. 



present article is to explore the role of a mutual interplay 
of quantum fluctuations and lattice vibrations in deter- 
mining the ground-state spin ordering, as well as, to ex- 
amine the influence of lattice vibrations on the magnetic 
behaviour in a close vicinity of critical points. Note fur- 
thermore that the amplitude of lattice vibrations in the 
most spontaneously long-range ordered magnetic materi- 
als is usually small enough below the critical temperature 
and hence, our further findings will be mainly confronted 
with the relevant results obtained for the rigid limit of the 
model under consideration [39] . 

Even though all the results derived in the foregoing 
section hold independently of the lattice geometry and 
also irrespective of whether exchange interactions are as- 
sumed ferromagnetic or antiferromagnetic, our particular 
attention will be hereafter restricted only to one represen- 
tative example of the spin-1/2 Ising-Heisenberg model on 
a doubly decorated square lattice with both ferromagnetic 
exchange constants J > and Ji > 0. Under this assump- 
tion, it is also reasonable to consider positive values of 
both magnetoelastic constants A > and Ai > 0, which 
ensure an enhancement (reduction) of the ferromagnetic 
interaction under the contraction (elongation) of the rel- 
evant bond length. For simplicity, let us also suppose the 
equality of the interaction parameters of the same charac- 
ter, i.e. J = Ji, A = Al, K = Ki, which will consequently 
reduce the total number of free parameters. Finally, all the 
interaction constants will be normalized hereafter with re- 
spect to a strength of the exchange interaction J so that 
the ratios k-QT/ J , A/ J and K/ J will mark a relative mag- 
nitude of the temperature, the magnetoelastic constant 
and the spring stiffness constant, respectively. 

First, let us take a closer look at the ground-state be- 
haviour. The typical ground-state phase diagram in the 
J/K—A plane is depicted in Fig. [21 Two first-order phase 
transition lines divide the ground-state phase diagram into 
three different regions corresponding to the classical fer- 
romagnetic phase (CFP), the quantum antiferromagnetic 
phase (QAP) and the third phase with immense ditor- 



tions of lattice bonds (FLP), which is beyond the scope 
of present approach (see for details the next two para- 
graphs). As one would expect, the perfect ferromagnetic 
alignment of all Ising and Heisenberg spins can be found 
in CFP, which however represents the ground state only if 
the exchange anisotropy A does not exceed one of its two 
critical values Ad or Ac2 ■ It is noteworthy that the highest 
critical value of A^{^^ = -\/3 corresponds to the rigid limit 
K/ J — > oo and one observes a very weak gradual decline of 
Aci with decreasing a relative strength of the spring stiff- 
ness constant K/ J that is only hardly discernible within 
the scale displayed in Fig. [21 On the other hand, the quite 
unusual QAP becomes the ground state above the critical 
value Aci for strong enough spring stiffness constants. It 
is worthy to recall that QAP surprisingly appears in spite 
of ferromagnetic character of all considered exchange cou- 
plings J = Ji > due to a mutual competition between 
the easy-plane XXZ Heisenberg interaction and the easy- 
axis Ising interaction (for more details see Ref. [39 ). The 
spin ordering inherent to QAP can be characterized by a 
perfect antiferromagnetic (Neel) long-range order of the 
Ising spins from the nodal sites of a doubly decorated lat- 
tice, which is accompanied with the symmetric quantum 
superposition of two antiferromagnetic spin states of each 
couple of the quantum Heisenberg spins from the same 
bond of a doubly decorated lattice [39 

It should be pointed out that CFP and QAP are the 
only two possible ground states, which are separated one 
from each other by the critical value A'^^^^ — a/3 when 
considering the rigid limit K/J^ oo [39 . For any fi- 
nite value of the spring stiffness constant, however, the 
investigated model enters FLP whenever the exchange 
anisotropy exceeds the second critical value Ac2- This lat- 
ter critical value Ac2 exhibits a rather steep decrease by 
shrinking the lattice rigidity, i.e. by lowering a strength 
of the spring stiffness constant K/ J . Consequently, QAP 
completely vanishes from the ground-state phase diagram 
for less stiff lattices with a relatively small spring stiffness 
constant, where the direct transition from CFP to FLP 
takes place since the critical value Z\c2 becomes smaller 
than Aci in this parameter space. 

The description of spin ordering inherent to FLP is 
however beyond the scope of present study, because one al- 
ways observes in FLP a rather sizeable mean displacement 
of the Heisenbergs atoms that is in contradiction with the 
harmonic approximation demanding small lattice defor- 
mations only. In order to provide a deeper insight into a 
magnitude of those lattice deformations, Fig.[31shows typ- 
ical zero-temperature changes of the average displacement 
of the Heisenberg atoms with the exchange anisotropy A 
for several values of the spring stiffness constant K/J. 
It indeed turns out that the mean displacement exhibits 
abrupt jumps if crossing the first-order phase transition 



^ Notice that the ground state of the model with the an- 
tiferromagnetic Heisenberg interaction J < is unique and 
it strongly resembles the spin ordering inherent to QAP with 
exception that the Heisenberg spin pairs reside the antisym- 
metric quantum superposition of two antiferromagnetic spin 
states rather than the symmetric one. 
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Fig. 3. Zero-temperature variations of the mean displace- 
ment of the Heisenberg atoms as a function of the ex- 
change anisotropy A for Ji/J ~ 1.0, A/ J — 1.0 and sev- 
eral values of the elastic constant K/J. 



lines from the ground-state phase diagram shown in Fig. [2] 
and it generally increases with the exchange anisotropy 
when keeping the same strength of the spring stiffness 
constant. The formula (j35p allows one to derive the fol- 
lowing ground-state values for the mean displacement in 
the individual phases 



(36) 



, , 2JAA^ - J,A, - A^Jl + {JAf 

^"''^"""^ AiK, + 2KUJl + iJA)^ ' ^''^ 

^^^^ (38) 



It can be easily understood from Eq. (|36|) that the average 
displacement in CFP equals zero just if the magnetoelas- 
tic coupling constants A and Ai are equal, otherwise the 
spontaneous distortion of the Heisenberg and Ising bonds 
takes place in CFP as well, li A > Ai, then, the con- 
traction of the Heisenberg bonds accompanied with the 
elongation of the Ising bonds occurs due to a larger mag- 
netoelastic energy gain of the former bonds, while the re- 
verse is the case if Ai > A. On the other hand, the mean 
displacement is always positive in QAP, which indicates 
shortening of the Heisenberg bonds at the expense of elon- 
gation of the Ising bonds on behalf of a more dominant 
magnetoelastic energy gain of the former bonds. The same 
holds true also for FLP, however, the mean displacement 
always gains in FLP rather high values, which are beyond 
the validity of our approach. 

Let us close our discussion of the ground state by in- 
vestigating the zero-temperature variations of the uni- 
form and staggered spontaneous magnetizations in de- 
pendence on the exchange anisotropy A for two different 
strengths of the spring stiffness constant K/ J ds shown in 
Fig. |4l It is noteworthy that the zero-temperature vari- 
ations of the uniform and staggered spontaneous mag- 
netizations depicted in Fig. |H[a) for the spring stiffness 
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Fig. 4. Zero-temperature variations of the sublattice mag- 
netizations (m^, ms) and the sublattice staggered magne- 
tizations {sA^ sb) as a function of the exchange anisotropy 
A iov Ji/ J — 1.0, A /J — 1.0 and two different values of 
spring stiffness constant: (a) K/J = 5.0; (b) K/J ~ 2.0. 



constant K/J ~ 5.0 are almost identical with the corre- 
sponding results of the perfectly rigid model {K/ J 00) 
given by 



CFP: [toa,tob,sa,ss] = 
QAP: [mA,mB,SA,SB]^ 



0,0, 



1 



Ji 



2'2^jf + {jAy 



(39) 
.(40) 



While the uniform sublattice magnetizations niA and uib 
of both Ising as well as Heisenberg spins achieve their sat- 
uration value in CFP and imply thus the classical char- 
acter of this spontaneously long-range ordered phase, the 
staggered sublattice magnetization sb of the Heisenberg 
spins is subject to the quantum reduction of the magneti- 
zation in QAP and only the staggered sublattice magne- 
tization SA of the Ising spins attains its saturation value. 
This latter finding evidently confirms the quantum nature 
of QAP, which is basically influenced by pronounced lo- 
cal quantum fluctuations that also suppress the staggered 
sublattice magnetization sb of the Heisenberg spins. Al- 
though the zero-temperature behaviour of the uniform and 
staggered magnetizations remains qualitatively the same 
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Fig. 5. The dependence of critical temperature on the 
exchange anisotropy A for Ji/J — 1.0, A/ J = 1.0 and 
several values of the spring stiffness constant K/ J. 



even for less rigid lattices, one may find an important fea- 
ture which directly follows from a comparison of Fig. |H[a) 
and (b) displayed for a more and less rigid lattices, re- 
spectively. It turns out that the local fluctuations of lat- 
tice spacings in conjunction with the local quantum spin 
fluctuations tend to enhance the quantum reduction in 
the staggered sublattice magnetization sb of the Heisen- 
berg spins, which originally comes from the competition 
between the easy-plane Heisenberg and easy-axis Ising in- 
teractions. Hence, it follows that the suppression of the 
staggered sublattice magnetization is the greater, the 
smaller is the spring stiffness constant K/ J, or the higher 
is the magnetoelastic constant A/ J. 

Now, let us proceed to a discussion of the critical be- 
haviour by investigating the finite-temperature phase di- 
agram shown in Fig. [5] in the form of critical tempera- 
ture vs. the exchange anisotropy dependence. It should be 
noted here that all the displayed critical boundaries have 
been obtained by solving numerically the critical condition 
l|23p . The finite-temperature phase diagram generally form 
two wings of the same critical line, which meet together at 
the critical value Ad separating in the ground state CFP 
and QAP. The left wing of a critical line actually repre- 
sents the phase boundary of CFP, whereas the right wing 
of a critical line represents the phase boundary of QAP. 
It is quite evident from Fig. [5] that the reduction in the 
lattice stiffness (rigidity) generally causes the suppression 
in the critical temperature regardless of whether CFP or 
QAP constitutes the ground state. It should be neverthe- 
less stressed that there is a quite substantial difference 
when comparing the size of this reduction in CFP and 
QAP. For instance, the decrement in the lattice stiffness 
by the factor of 50 suppresses the critical temperature of 
CFP at most by 10%, while the same decrement in the lat- 
tice stiffness may cause in QAP much more substantial di- 
minishing of the critical temperature or even the complete 
loss of spontaneous long-range ordering (e.g. compare the 
curves K/J = 100 and 2 displayed in Fig. [5]). The attenu- 
ation in the lattice rigidity may be thus responsible for an 
absence of the spontaneous ordering inherent to QAP (i.e. 
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Fig. 6. Temperature dependences of both spontaneous 
sublattice magnetizations (to^, ms) for Ji/J = 1.0, 
A/ J = 1.0, two different values of the exchange anisotropy 
A — 1.0, 1.5 and two different values of the spring stiffness 
constant: (a) K/J = 5.0; (b) K/J = 2.0. 



the right wing of a critical line) as it is demonstrated in 
Fig.[S]on the particular case with K/J = 0.5. This obser- 
vation would suggest that the lattice distortions basically 
reinforce the effect of quantum fluctuations in destroying 
the spontaneous long-range order of the quantum nature 
(such as QAP), while the disappearance of classical spon- 
taneous long-range order (such as CFP) is substantially 
less affected by a presence of the lattice deformations. 

Last but not least, we will turn our attention to the 
temperature dependences of the order parameter, which is 
the uniform spontaneous magnetization for CFP and the 
staggered spontaneous magnetization for QAP. Fig. [5] il- 
lustrates typical thermal variations of both uniform spon- 
taneous sublattice magnetizations for two different values 
of the exchange anisotropy and two different values of the 
spring stiffness constant. It is worthy to notice that the 
sublattice magnetization niA of the Ising spins is always 
more robust with respect to the thermal fluctuations than 
the sublattice magnetization ms of the Heisenberg spins 
even if both spontaneous sublattice magnetizations disap- 
pear in a vicinity of the critical temperature with the same 
critical exponent /3e = 1/8 from the standard Ising univer- 
sality class. It can be also clearly seen from the compari- 
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son of Fig. [Blja) and (b) that the reduction in the lattice 
stiffness does not change basically the temperature de- 
pendences of the spontaneous magnetizations, but it only 
causes a small gradual shift of the critical temperature 
towards lower values. 

The relevant temperature behaviour of both staggered 
spontaneous sublattice magnetizations in QAP is illus- 
trated in Fig. [71 While the critical exponent determin- 
ing the variation of both staggered sublattice magnetiza- 
tions in a close vicinity of the critical temperature remains 
the same Pe — 1/8, the drop in the lattice rigidity has 
much more obvious effect upon the thermal dependences 
of the staggered sublattice magnetizations. It can be in- 
deed readily understood from the comparison of Fig. [71[a) 
and (b) that the staggered sublattice magnetization sb 
starts from much lower initial value when decreasing the 
lattice stiffness, because it becomes subject of a stronger 
quantum reduction of the magnetization due to the inter- 
play between the local quantum spin fluctuations and the 
local fluctuations of lattice spacings. In addition, one also 
observes a more substantial reduction in the critical tem- 
perature by lowering the spring stiffness constant, which 
comes from a mutual interplay between the temperature- 



induced spin fluctuations, the quantum spin fluctuations 
as well as the distortion of lattice bonds. 



4 Conclusion 

The present article deals with the magnetoelastic prop- 
erties of the spin-1/2 Ising-Heisenberg model defined on 
doubly decorated planar lattices whose decorating atoms 
may oscillate around their equilibrium lattice positions af- 
ter taking into account the elastic energy penalty of each 
distorted bond within the harmonic approximation. It has 
been shown that after tracing out the spin degrees of free- 
dom of the Heisenberg spins and integrating out the lat- 
tice degrees of freedom it is possible to establish a sim- 
ple mapping equivalence between the investigated model 
system and its corresponding spin-1/2 Ising model on un- 
decorated rigid lattice using the generalized decoration- 
iteration transformation [52153154155] . With the help of 
this mapping method, we have examined the ground-state 
and finite-temperature phase diagrams, the magnitude of 
the mean displacement, as well as, the temperature de- 
pendences of the order parameters. 

To shed light on the influence of the lattice vibrations, 
our particular attention has been devoted to the most 
significant differences between the magnetic behaviour of 
the spin-1/2 Ising-Heisenberg model on the perfectly rigid 
doubly decorated lattice [39] and respectively, the dou- 
bly decorated lattice whose decorating atoms are relaxed 
from the condition of a perfect rigidity. It has been demon- 
strated that the mutual interplay between the quantum 
spin fiuctuations and the lattice deformations enhances 
typical quantum features like the quantum reduction of 
the magnetization in the ground state of quantum phases 
such as QAP, while it does not affect the respective ground- 
state behaviour of classical phases such as CFP by no 
means. As far as the magnetic behaviour at finite tempera- 
tures is concerned, it also turns out the fluctuations in lat- 
tice spacings act more efficiently in conjunction with ther- 
mal and quantum spin fluctuations in destroying a spon- 
taneous long-range order of the quantum phases, whereas 
the influence of the lattice distortions in promoting ther- 
mal and quantum fluctuations by disturbing a sponta- 
neous long-range order of the classical phases is much less 
efficient and conspicuous. 

Finally, let us briefly mention that theoretical treat- 
ment of lattice vibrations within the framework of clas- 
sical harmonic oscillators rather than quantum harmonic 
oscillators represents perhaps the most important draw- 
back of the proposed model. It should be nevertheless 
mentioned that the approach presented in this paper can 
be rather straightforwardly generalized to the spin-1/2 
Ising-Heisenberg model on doubly decorated planar lat- 
tices whose lattice vibrations will be treated as quantum 
harmonic oscillators as well. Our preliminary study of 
this latter more sophisticated model reveal the completely 
same results for the ground-state and finite-temperature 
phase diagrams, as well as, the relevant behaviour of the 
order parameter due to the same form of the magnetic 
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part of the bond partition function |72j . The most fun- 
damental differences between the Ising-Heisenberg model 
with the lattice vibrations described at the classical and 
quantum level, respectively, thus appears in a respective 
behaviour of basic thermodynamic quantities (such as in- 
ternal energy, entropy or specific heat) that are also af- 
fected by the elastic part of the bond partition function. 
For instance, the former model apparently exhibits arti- 
ficial temperature dependences of the specific heat owing 
to the anomalous lattice contribution of the heat capacity, 
but we have not dwell on investigating this aspect of our 
model in the present paper. 

The financial support of this work provided by tlie research 
foundation FAPEMIG under the grant number BPV-00088-10 
is gratefully acknowledged. 
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